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ABSTRACT 

In this paper we study the angular momentum properties of simulated dark matter 
halos at high redshift that likely host the first stars in the Universe. Calculating the 
spin distributions of these 10^ — 10^ halos in redshift slices from z = 15 — 6, we 
find that they are well fit by a log-normal distribution as is found for lower redshift 
and more massive halos in earlier work. We find that both the mean value of the 
spin and dispersion are largely unchanged with redshift for all halos. Our key result 
is that subsamples of low and high spin 10^ and 10^ halos show difference in 
clustering strength. In both mass bins, higher spin halos are more strongly clustered 
in concordance with a tidal torquing picture for the growth of angular momentum in 
dark matter halos in the CDM paradigm. 
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redshift ~ galaxies: formation 
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1 INTRODUCTION 

In the standard ACDM paradigm of structure formation, 
dark matter over-densities accrete mass, forming deep po- 
tential wells into which baryons condense and cool to even- 
tually form stars. The details of the cooling and star 
formation process are quite complex, and the physics of 
the first episode of star formation and that of subse- 
quent generations is not completely understood. There have 
been detailed numerical studies of the formation of the 
first stars which suggest that th e masses of these stars 
were skewed to the massive end jAbel. Bryan, fc NormarJ 



2002l:lBromm. Coppi. fc Larsonl2002l : IOmukai fc Pallall2003l : 



O'Shea fc NormanI 120071 ) . Numerical simulations are chal- 
lenging owing to the large dynamical range needed as well 
as the complexities of hydrodynamical modeling required 
to capture the essential physics. Current models of galaxy 
formation in the early Universe are able to make some 
simplifying assumptions to make the si mulations of sub- 
sequent star formation easier (s e e, e.g., Grief et al] 20081: 



Johnson. Grief, fc Bromm 2p08|: Ricotti. Gnedin. fc Shulll 



2002al lbl: IWise fc Abell l2007l : IWise. Turk, fc Abell |2008| ). 

However, much work still needs to be done to fully uncover 
the details of galaxy formation, including a better under- 
standing of the properties of the dark matter halos and how 
the dark matter properties affect baryonic structure forma- 
tion. 

In numerical simulations, three basic properties of dark 
matter are followed; mass, position, and velocity. From these 



three properties, much has been learned about the buildup 
and assembly of early dark matter halos, including halo 
mass f unctions, galaxy c lustering, halo shapes, and halo 
spins dRe ed ct al. 2007; Jang-Condell fc Hernquist '200l|; 
Reed et al. 2003, 2 008 : Maccio, Button, fc van den Bos^ 
2008h . Refinement in our understanding of dark matter halos 



is needed to gain insight into how the baryonic component 
of galaxies couples with the dark matter (e.g., angular mo- 
mentum transfer, stellar and AGN feedback) . A key missing 
piece is the role of angular momentum in the formation and 
modulation of structure formation. This paper focuses on 
the angular momentum and clustering properties of dark 
matter halos in the early {z — 15 — 6) Univers, roughly 
2 Gyr after the Big Bang in the concordance model. The 
goal of this work is to understand more about the early, low 
mass dark matter halos that likely host the earliest stars and 
galaxies in a statistically large sample. 



2 ORIGIN OF ANGULAR MOMENTUM IN 
DARK MATTER HALOS 

In the standard ACDM paradigm, angular momentum is ac- 
quired by dark matter halos due to tidal torquing by their 
neighbors l|Hovlelll949l ). During the linear phases of proto- 
galactic evolution the angul ar momentum grows to first or- 
der in proportion to time t. iPeebled (|l969l ) found from the 
study of collapsing spherical regions that the growth in the 
angular momentum J occurs in the second order only as 
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t^^^. In such regions growth occurs purely as a result of con- 
vective effects on the bounding surface. This was calculated 
using linear theory to find the growth rate of the angular 
momentum within a co-moving sphe rical region of an ex - 
panding Einstein de-Sitter Universe. iDoroshkevichI (|l970l ) 
showed that while the angular momentum of a proto-galaxy 
did grow cx t for a flat Universe, the Peebles result was 
a consequence of the i mposed spherical symmetry and was 
verified numerically bv lWhitj (il984i ) . We briefiy outline the 
process of angular momentum acquisition by a dark matter 
halo in what follows. 

Consider the growth of fiuctuations in an expanding 
Friedmann Universe filled with pressureless matter with den- 
sity, p{f,t), and let 



5{r,t) 



-1; S{x,t) =b{t)5o{x). 



(1) 



po(f,t) 

The local over-density can be written as a separable function 
of time and a co-moving coordinate x related to r by r = 
a{t)x, where a{t) is the cosmological expansion factor. The 
trajectory of a dust particle in linear theory is given by: 



riq, t) = a{t) x{q, t) = a{t) [q - b{t) V<l>{q)] 



(2) 



where q is the Lagrangian coordinate defined as the x posi- 
tion of the particle as t ^ and (p the peculiar gravitational 
potential. If (^6^'^ ^ 1, a condition that is satisfied over the 
period of time preceding the collapse only when the initial 
density field has a coherence length of the proto-galactic 
scale. In general, this will cease to describe the detailed evo- 
lution of the matter distribution long before galaxies form. 
The spin angular momentum of the material that makes up 
a proto-galaxy can be written as. 



J{t) — poa' / {x — x)Xxd''q 
Iv, 



(3) 



Since the leading term in x is first order and is parallel to 
x, the 2nd order term gives the non-zero contribution to the 
above integral. Writing it explicitly as an integral over the 
surface (using Gauss's theorem), we have: 



(4) 



where T is the tidal tensor and describes the local deforma- 
tion at q and / is the inertia tensor of the matter in Vl- 
This tensor product is used to calculate the tidal torque on 
an extended body in the tidal field. Angular momentum is 
acquired as the first order tidal field couples to the zeroth 
order quadrupole moment of the irregular boundary of the 
proto-galaxy. The principal axes in general do not coincide 
as T defines the shape and disposition of the neighboring 
perturbations and / depends only on the shape of the proto- 
galaxy. This outlines in principle the origin and offers insight 
into the calculation of angular momentum of simulated dark 
matter halos. 



3 METHOD 

We ran a series of N-body simulations to follow the growth 
of dark matter halos from z ~ 100 down to z = 6. In order 
to study halos that likely host Pop III stars, we chose the 
particle mass such that a 10® M© dark matter halo would 
have ~ 100 particles. For 512"^ particles, this requirement 



sets the box size at 2.46Mpc/h and the particle mass at 
Mdm = 7.3 X 10"^ M0. The initial conditio ns were gen- 
erated using a parallelized version of Grafic IjPrunet et al.l 
,2008.), which calculates the Gaussian rand om field for the 
dark matter particles. We used Gadget-2 ISpringell [iooH ) 
to follow dark matter particles down to a redshift of 2; = 6, 
with output s napshots at z = 15, 1 2, 11, 10, 9, 8, 7, 6. We used 
the WMAP3 (Sper gel et alll2007l ) cosmological parameters 
for the runs presented here, [Qm ,nA,^b, h,n, as] ~ [0.238, 
0.762, 0.0416, 0.732, 0.958, 0.761]. 

For each snapshot, we used a friends-of- friends code to 
find individual halos, using a linking length of 0.2 times the 
initial particle spacing. Once the halos were identified, each 
particle in the halo was teste d to see if it was actually bound 
to the halo. We used SKID (IStadelll200ll ') to do the unbind- 
ing. SKID finds the potential and kinetic energies for all 
particles in a given halo and removes the most unbound par- 
ticle from the halo. Successive iterations are performed until 
all particles are either bound or there are no more particles 
in the halo. Without unbinding, there was excess power in 
both tails of the spin distribution. 

Once unbound particles were removed, for each halo 
that had more than 50 dark matter particles, we calculated 
the dimensionless spin parameter, 



A = 



J\E 



1/2 



GM5/2 ' 

where J is the total angular momentum, E the total en- 
ergy (kinetic plus potential), and M the total mass of the 
halo. For the most massive halos, with more than 10^ parti- 
cles, we used the hier archical force calculator Treecode vl.4 
dBarnes fc Hutlll986l ) to get the potential energy, inducing 
an error of 0.2 % in the potential energy as compared to 
direct summation. 



4 RESULTS 

Figure [T] shows the spin distribution at four redshifts 
{z = 6, 10, 12, 15). A s has been found for more massive 
halos at low redshif t (iBett et al. 20071: Bailin fc Steinmet3 
20051: ICole fc Laceyl Il996l : ISteinmetz fc BartelmannI Il995r 



Warren et all 1 19921 e.g.,), the distribution is log-normal, 
with a slight excess at l ow sp in values. At z = 10, 
IJang-Condell fc Hernguistl (|200ll ) also found a log-normal 
distribution, with (A) = 0.033 and a = 0.52, while we 
find (A) = 0.04 and a = 0.55 at the same redshift. These 
differences may be due to the increased resolution, vol- 
ume, and number of particles i n our study. In addition, 
Ijang-Condell fc Hernguistl (|200ll ) used a different halo find- 
ing algorithm (SKID's density maxima finder), and found a 
Press-Schect er mass fun c tion. However, our mass function 
matches the iReed et al.l (j2007? ) modification to the Sheth- 
Tormen mass function. Finally, the differences between their 
results and our work may be due to different cosmological 
parameters chosen, particularly ag. In future work, we plan 
to study the effect of mass resolution on the derived angular 
momentum properties as well as the effect of changing as- 

The right column in Figure [1] shows the trend in the 
mean and standard deviation of the fitted log-normal over 
redshift. The error bars show the la estimated errors on the 
mean and standard deviation of the fitted log-normal. We 
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Figure 1. Spin parameter distribution for all halos at four redshifts. The solid curve is a log-normal fit to the data. The right two panels 
show the change in mean A and the standard deviation for the log-normal fit for all redshift outputs used in this study. 



found only slight variation in the mean spin, with A = 0.0378 
at z = 15 increasing to A = 0.0426 at 2 = 6. The small 
increase in spin over time is most likely due to stronger tidal 
torques from more massive halos that collapse at later times. 
The standard deviation decreases slightly from cr = 0.573 at 
z = 15 to a — 0.553 at 2 = 6. 

We show in Figure [2] the spin distribution for halos with 
mass within ±0.2 dex of lO'^ M0 and lO'' M0 at the same 
redshifts as Figure [1] For the 10^ Mq halos, we find that (A) 
does not monotonically increase with redshift as it did for 
the whole sample, and (A) is systematically larger than the 
sample mean at any given redshift. There is also a sharp 
increase in (A) at 2 = 15, possibly because these halos were 
still in the process of forming and may have significant in- 
falling mass clumps. However, there are only 55 halos in 
this mass bin at z = 15, making it difficult to reliably fit 
the sample with a log-normal curve. The 10^ Mq halos show 
spin distributions more similar to the entire sample, with a 
slowly increasing mean spin from 2 = 15 to 2 = 6. 



5 CORRELATION FUNCTION OF 

HIGH-REDSHIFT DARK MATTER HALOS 

We also studied the clustering properties of dark mat- 
ter halos selected by their angular momentum. It is 
well known that more massive halos are strongly clus- 



tered 



(e^ 



Mo fc Whit3 I2OO2I; [Sheth fc TormenI 



iLemson fc Kauffmann|[l999h . Lemson fc K aufi^mann (199^) 
found no correlation between A and environment for 



;ly clus- 
n! Il999|; 



halos with Mtot > 10 ^ ^ Mr at 2 = 0, (see also 
Cervante s-Sodi et al.ll2008l : iMaccio et al]l2007l : iBerta et all 



Bct t et alJbOOTl '). This diagnostic is potentially im- 



2008; 

portant as a halo of mass M with larger angular momentum 
may take longer to collapse and therefore longer to form 
stars than a lower spin halo of the same mass. Thus, low 
spin halos could preferentially host star formation earlier. 
This in turn may allow feedback processes such as ionizing 
photons, stellar winds, and supernovae to become important 
earlier than in their high spin counterparts. These are im- 
portant consequences for galaxy formation, and the detailed 
role of the spin of the halo on the collapse of baryons will 
be explored in a subsequent paper. 

We calculate the two-point correlation function, ^(r), 
using the standard definition. 



DD{r)RR{r) 
i?D(r)2 



- 1, 



(5) 



where DD{r) is the number of halo-halo pairs within a sepa- 
ration of r — dr/2 and r + dr/2, RR{r) the same for a random 
distribution, and RD{r) is the cross-correlation between the 
data set and the random distribution. The quantity ^(r) rep- 
resents the excess probability of finding a halo at distance r 
compared to a random distribution. 

Figure 13] shows the correlation function, £,(r), for halos 
grouped by spin parameter for a given mass. The top panel 
is for halos with mass 10^ Mq, and the bottom panel for 
10^ Mq halos. We separated the halos out by mass to delin- 
eate the effect of spin from that of mass. We first computed 
the overall correlation function for all halos in the given mass 
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Figure 2. Spin parameter distribution for two halo mass bins, 10^ Mq (crosses) and 10^ Mq (diamonds). The solid (10^ Mq) and dashed 
(10^ Mq) curves are log-normal fits to the data, and the right two panels show the change in mean A and the standard deviation for the 
log-normal fit for each of the mass bins. 



range, and then used two sub-samples based on their spin 
parameter. The cuts in spin parameter were chosen so that 
one third of the halos were in the high and low spin bins 
respectively. The error bars show Poisson errors in DD{r), 



RR{r) 



RD{r) 



When counting halo pairs for the spin cut samples, we em- 
ployed two methods. First, we counted only pairs where both 
halos were in the given bin, and these results are shown in 
Figure [3] We also counted the pairs where only one halo 
was in the high or low spin bin. This method of counting 
produced similar results as the first method. 

We find that there is a distinction in the correlation 
function when the data are separated by the spin param- 
eter. High spin halos are more strongly clustered than low 
spin halos, similar to the result found bv lBett et al.l (|2007l ) 
in the Millennium Simulation. We also find that the differ- 
ence in clustering is preserved over a range of redshifts. The 
excess clustering is due to the fact that in a denser environ- 
ment, halos feel stronger tidal torques, and thus have larger 
angular momentum and spin parameters. The correlation 
function can be fit to a power-law, ^(r) — (r/R)'^, where R 
is the correlation length. We find the high spin halos have a 
correlation length on average 25% larger than the low spin 
halos over the various mass and redshift slices used here. 
The slope, 7, did not show any correlation with spin, mass, 
or redshift. 



6 DISCUSSION 

In this paper, we have examined the angular momentum 
distribution of simulated high redshift dark matter halos 
and their correlation functions. While the angular momen- 
tum distributions for this population are well fit by log- 
normal distributions as is the case for lower redshift and 
more massive halos, we find that the clustering properties 
of the high redshift population depend on the value of spin. 
For a given mass bin selecting by spin, we find that the cor- 
relation function is higher for high spin halos. This is an 
important trend and appears to be robust. Although in this 
study we have restricted ourselves to the analysis of dark 
matter only simulations, the different correlation lengths of 
high and low spin halos are likely to have an important im- 
pact in feedback from these early galaxies. Pop III stars 
appear to h ave a large impact on their environment due 
to ra diative (| Johnson. Grief, fc Brornrn 2007: Whalen et al.l 
l2008h and supernova (|Grief et al...2007. : .Whalen et al.ll2008l ) 
feedback. 

W e note here that in a recent paper [O'Shea fc Normal] 
(|2007l ) tracked the spin parameters of halos that formed the 
first stars in 12 cosmological random realizations and did not 
find a correlation between halo spin and collapse time. The 
resolution and analysis pursued by them differ from ours as 
they zoom in and track (by explicit selection) the properties 
of the dark matter halo that hosts the first star. Our analy- 
sis is distinct by construction as we are measuring the spin 
parameters of an ensemble of dark matter halos and arguing 
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Figure 3. Correlation function for 10^ M© (top panel) and 10^ M© (bottom panel) halos as a function of comoving separation, r. Each 
plot shows the correlation function for all halos in the mass bin (solid), as well as for a cut based on the halo spin parameter. The cuts 
were chosen so that one third of the halos would fall in the high (dotted) and low (dashed) bins. 
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that the higher spin halos might in general commence the 
collapse and formation of its first star later. As shown in 
Figure [21 the distribution of halo spin for the high mass bin 
at 2 = 15 is shifted to higher spins when compared to the 
general population. Thus, the larger halos that would host 
Pop III stars will h a ve lar ger spins, but the largest spin that 
lO'Shea fc NormanI (|2007l ') followed was for A « 0.1, which 
is at the half- max of the log-normal fit. One halo with this 
spin may not be representative of all high spin, massive ha- 
los, which may show a delayed baryonic collapse. While this 
is beyond the scope of this paper, we intend to pursue the 

consequences of our results for baryonic collapse. 

Our results are in agreement with what iBett et aP 

found for the correlation of dark matter halos at 
2 = in the Millenium Run. lBerta et al] ()2008h found trends 
in the SDSS spectroscopic sample of local galaxies relating 
the estimated dark matter spin parameter with the galaxy's 
baryonic properties, such as stellar mass and recent star for- 
mation histories. Thus, angular momentum may be a key to 
determining the epoch of formation, size, and structure of 
the stellar and gaseous disks in galaxies. In the early Uni- 
verse however, it remains unclear how much effect angular 
momentum will have on the first generation of galaxy for- 
mation, although our results are suggestive. 
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